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Abstract 



^ . KdV6 equation can be described as the Kupershmidt deformation of the 

O I KdV equation (see ^005, Phys. Lett. A 372: 263). In this paper, starting from 

the bi-Hamiltonian structure of the discrete integrable system, we propose a 
cn ■ generahzed Kupershmidt deformation to construct new discrete integrable 

'!::j- \ systems. Toda hierarchy, Kac-van Moerbeke hierarchy and Ablowitz-Ladik 

C3 ' hierarchy are considered. The Lax representations for these new deformed 

systems are presented. The generalized Kupershmidt deformation for the 
discrete integrable systems provides a new way to construct new discrete in- 
tegrable systems. 
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Key words: Kupershmidt deformation; bi-Hamiltonian systems; discrete integrable 
systems. 



1 Introduction 

In recent years, the deformations of integrable systems attracted lots of attention 
[T]-[I7]. It is known that one can construct a new integrable system from a bi- 
Hamiltonian system [21 [3]. KdV6 equation is derived by A. Karasu-Kalkani, et al.. 
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To] by means of the Painleve analysis. KdV6 equation is written as 

Ut = Uxxx + Quux - Ux, (1.1a) 

^xxx + 4:UUx + 2uxUJ = 0, (1-lb) 



which could also be seen as the nonholonomic deformation of KdV equation |10] . 
Many authors studied the integrability properties of KdV6, such as zero-curvature 
representation, bi-Hamiltonian structure, conserved quantities, multisolitons and so 
on [ID]- [15]. Kupershimidt found that ( II. ip can be converted into 

Ut = A^) - J{uj), (1.2a) 

ou 

K{uj) = 0, (1.2b) 

where 

j = d = Q^^ K = d^ + 2{ud + du), 



„2 



are two standard Hamiltonian operators of the KdV hierarchy and H^ = u — f is 
the Hamiltonian function. In general, for a bi-Hamiltonian system 



ut^ = JG^^^) = K{- 



6u 6u 

there is a nonholonomic deformation of bi-Hamiltonian system [1] 

ut^ = J{^^) - Aoo), (1.3a) 

ou 

K{u) = 0, (1.3b) 

which is named as Kupershmidt deformation of bi-Hamiltonian systems. 

In [T5|, we proved that the Kupershmidt deformation of the KdV equation could 
be seen as the Rosochatius deformation of the KdV equation with self-consistent 
sources, and its bi-Hamiltonian structure are obtained. We also proposed a gener- 
alized Kupershmidt deformation (GKD) for some bi-Hamiltonian system with pure 
differential Hamiltonian operators in [16] and [T7]. For an integrable system, we 
know that its two Hamiltonian operators have the relation 

where -j^ can be taken from the spectral problem of the system and 7^ = Xj or 
7j = A^ depending on which integrable system we are studying. For the integrable 
system with self-consistent sources, we have eigenfunctions and adjoint eigenfunc- 
tions as components of the sources and they will satisfy the spectral problem and 
the adjoint spectral problem corresponding to Xj. In the GKD case, we assume that 
the nonholonomic term satisfy {K — ljJ){~si) ~ ^ instead of the spectral problem 
so that they will have more freedom as they satisfy higher order equations. For the 



continuous case, we have found that the GKD of the integrable system will have the 
form [Ml [IT] 

(K - 7,^)(^) = 0, J = 1, 2, . . . , AT, (1.4b) 

where •jj equals to Xj or Xj depending on the integrable system. 

The integrability of the generalized system is proved by constructing its Lax 
representation. The present paper focuses on the GKD for the discrete integrable 
systems. We take Toda lattice hierarchy, Kac-van Moerbeke hierarchy and Ablowitz- 
Ladik hierarchy as examples to construct the generalized Kupershmidt deformation 
for these systems and their Lax representation are obtained to show that the new 
systems are integrable generalization of the original discrete systems. 

Zhou introduced a concept of mixed hierarchy of soliton equations combining the 
negative order equations to the positive order equations and showed that the Kuper- 
shmidt deformed systems are just special members in the mixed soliton hierarchies 
|14] . The mixed hierarchy of soliton equations [H], which has one nonholonomic 
term, are some special nonholonomic deformed equations. Our GKD admit A^ non- 
holonomic terms. 

The paper is organized as follows. In section 2, we construct the GKD for the 
Toda hierarchy. Section 3 treats the GKD for the Kac-van Moerbeke hierarchy. 
In Section 4, we will obtain the GKD for the Ablowitz-Ladik hierarchy. Some 
conclusions will be given in Section 5. 



2 The generalized Kupershmidt deformation of 
Toda hierarchy 

Let us recall some fundamental concepts of the discrete integrable system. Assume 
/ = f{n,t) for ra G Z and t G M. The shift operator E and difference operator D 
are defined as 

{Ef){n) = f{n + 1), {E~'f){n) = f{n - 1), 

iDf){n) = {E-l)f{n), n e Z, (2.1) 

and /(j) = E^f, j e Z. 

Consider the following discrete eigenvalue problem for the Toda hierarchy [18] 

^(-1) = f/(u,t;,A)^, (2.2) 

where 



First, we consider the following stationary zero-curvature equation for generating 
function F, 

T^~^^U-UT = 0, (2.3) 

where 



+00 +00 



tti hi 



and Oj, 6j, Cj are functions of n and t. In general, the recursion formula is 

where the initial data is chosen as oq = |, 60 = Cq = 0. Define the modification 
matrix 

Am := diag{bm+i + S,S), 

where 6 is an arbitrary constant. Let 

*t. = Kn^. (2.4) 

where 

Kn := (A'"F)+ + A„, = J2 r.A™-* + A^^. 

i=0 

Then the compatibility condition of (12. 2 p and (12. 4 p gives rise to the zero-curvature 
representation of the Toda lattice hierarchy 

Ut^ = UVm - Vi-'^U. (2.5) 

For m = 1, we have the famous Toda equation 

vt = v{u^^^^ -u), (2.6a) 

ut = v- v^^^ (2.6b) 

The bi-Hamiltonian structure of the toda hierarchy reads (see, e.g., [19]) 

where 

^(T) _ f -vEv + vE^^v -vu + vE^^u \ , , . 

~V vu-uEv vE-'-Ev J' ^ ' 

are two standard Hamiltonian operators. 

Now we derive -jf and -^ from (12. 2 p and its adjoint problem with A = Xj. 
Define L = v^^'E + u-\- E~^ and its adjoint operator is L* = vE~^ + u + E. Then 



we have Lipj = Xjipj, L*(j)j = Xj(f)j. We assume that L, Xj, ipj and (pj have a 
perturbation e. So we can assume that they are functions of e. Furthermore we can 
properly choose (pj to make < ipj,4>j >= 1, where the inner product is defined as 
< f,g >= J2n /('^)5'('^)- With all these assumptions we have 

^(^j(^)) = ^(< >'J^J^<Pj >) = ^(< LtpjAj >) 

dL dijjj ^ ^ , d(f)^ 

=< ^^i' 0i > + < L-^, 0J > + < Ltjjj, -^ > 



dv^^'^ (I) du , ^ dip 4 ^^ ^ ^ ^ d(j)j 



dv^^^ II) du , ^ . d f , , X 

dv^^^ , (1) du , 



So we have 



,1 = 1 M^"^ 






- ' ^ ^i(t). 



Take 7j = Aj, fll.4p gives rise to the following new generalized Toda lattice 
hierarchy 



\ / (-1) / ^ ^ 

^ \ _ j{T)l\ (Kn+ll^' 



u J, \ -h 






jiT,^ a ,,/v X_y i.,_^^ ^^ ^2.8a) 



(i^(^) - A, J(^)) (^ ^A; j=0, j = l,2,...,iV. (2.8b) 

It is noted that in our new deformed system, ipj and (pj are no longer the eigenfunc- 
tions and adjoint eigenf unctions of the spectral problem, but are solutions of fl2.8bp . 
After simplifying fl2.8bp . we can find that 



+ u(-i)v]"'^0j"'^ - A.^j-'^^J."^) + X.^^^. = 0, (2.9a) 

VU1pj(pj — UV^^'ijjj '(pj + V^~ (pj 

- v^'^^pf^f + X.v^'^^f^j - Xjvijj^'--'^ = 0. (2.9b) 

Denote fj = uipj + v'^^^ipj + ip~ — Xjipj and Qj = u(pj + v(pj + (pj Xj(pj, the 
above equations will lead to 

<^./. - ^^M"'^ = 0, (2.10a) 

^-1)/^. _ ^{i)^«^^, = 0, (2.10b) 



which gives that 

93 = =^^^7TT-- (2-llb) 

where /ij are some arbitrary constants. Therefore, we make the simphfication of 
(I2.8bp and find that our new constraints are the original spectral problem with 



some addition terms with some arbitrary constants /ij. the GKD of Toda lattice 
hierarchy is 



V 



uij, + v^^Hf + ^j"'^ = \i^j + _ ^'-,_iy (2- 12b) 

nh + ^^-'^ + '^f = hh + ,,,-^'jm. - (2-12C) 

We can find that in the GKD of the Toda hierarchy, the constraints is the spectral 
problem and its adjoint problem with some nonholonomic terms with arbitrary 
constants /ij. 

When 777, = 1, the GKD of Toda equation reads 



1—1 1—1 


(2.13a) 


TV N 




7-1 1-1 


(2.13b) 


„„/, 1 „,(1),/,(1) 1 „/,(-!) - \ ./. 1 ^i 


(2.13c) 


vi)j(t)j(t)) 




(2.13d) 


and its Lax representation is 




(w(^)E + m + E-^)^ = A^, 


(2.14a) 


N (1)7 


(2.14b) 



i=i 



The compatibility of the system (I2.14p under the condition (I2.13cp and (]2.13dp will 



give (12.13b ) and (12.13b ). It is an interesting question to find a Lax pair with /ij inside 
to give the system (I2.13p . This will be studied in future. When we take \j = 0, 
(I2.13P reduces to the Kupershmidt deformation of the Toda equation. When we 
take /ij = 0, (I2.13P reduces to the Toda equation with self-consistent sources [20] . 



3 The generalized Kupershmidt deformation of 
Kac-van Moerbeke hierarchy 



Consider the eigenvalue problem for the Kac-van Moerbeke hierarchy 

{E + vE-^)iJj = Xijj. (3.1) 

Its adjoint equation reads 

{E-^ + Ev)(l) = X(j). (3.2) 

We can also write the eigenvalue problem in matrix form 

To derive the hierarchy of evolution equations associated with the eigenvalue prob- 
lem, we frst solve the stationary discrete zero-curvature equation 

(3.4) 
(3.5) 

(3.6) 



r«[/ -UT = 0. 




Let 

r= " ^ =Vr.A- = y "^ ^^ A-^ 


i 


and take the initial data as 




ao =2' ^0 = 0, 6i = -1. 




We can find the recursion formula 




&2i = C2i = a2i+l =0, z = 0, 1, . . . 




C2J+1 — —V02i_^i, 




b2i+i = -(47^^ +a2i), 




cS+l + ^hi+l = a[%2 - «2i+2. 




Set the auxiliary linear problem as 




V't,™ = '^2mV^ 




with 




V2m = (rA2™)+ + A2„ 











(3.7) 



(3.8) 

Then the compatibility condition of (13. 3p and (13. 7p gives rise to the zero- curvature 
representation of the Kac-van Moerbeke hierarchy 

Ut..=villU-UV2^ (3.9) 



When m = 1, we have the Kac-van Moerbeke equation 

y^ = y^yi-^) ^ yWy (3.10) 

The Hamihonian operators of Kac-van Moerbeke hierarchy are defined as 

jiKvM) ^ ^(^-1 _ ^^^^ (3_^^g^) 

J^iKvM) ^ ^^^^-1 ^ ^(-1)^-1 + ^(-1)^-2 _ ^,^ _ ^(1)^2 _ ymE)y. (3.11b) 

The Kac-van Moerbeke hierarchy has the bi-Hamihonian form 

where if 2™ = — 1?^^ is the Hamiltonian function. It is not difficult to find that 



'2m 2m 



§ = ^]-'V.. (3.13) 

The new generahzed Kac-van Moerbeke hierarchy is constructed as follows 

AH- ^ 

v,^ = ji^"M)^J^ + J(^^^^) 5:(^j-^)0,), (3.14a) 

(^(i^'^M) _ a2j(^-^^))(^(-i)0^.) = 0. (3.14b) 



We assume that 



{E + vE-^)ijj = Aj^j + /,-, (3.15a) 

(E^^ + Ev)(l)j = \j(t)j + ^j. (3.15b) 

By simplifying fl3.14bp we have 

+ .(-^)(A,^(-^) - i,Y%Y'^ + .("^)(A,05-^) - v-^,)f^-'' 

+ ^^"'^^]"'Vi"'^ - ^^'^^fV, + (^^'^ - v^-^^){i,,^\r^^ - ^"'^0,) = 0. (3.16) 
To calculate (1 3 . 1 5 al) x (j) j — fl 3 . 1 5b p x ^ j . we will find that 

Therefore we make the assumption that 

^;0(.-i)/j-i) = v^P^g^, (3.17) 

which will lead to 

^j9j = ^Pj^^~'^-vtlj\-'^^,. (3.18) 
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Substituting f l3.18p into f l3.16p and with the assumption 

^;«0«/,=^^]"'^^i, (3.19) 

we will find (13.1 6p holds. With the assumption (I3.17p and fl3.19l) . we can solve for 
fj and Qj as 

fj = ^^^^TT> (3-20a) 

g, = !^ (3.20b) 

where Hj are some arbitrary constants. So the GKD of Kac-van Moerbeke hierarchy 
is 

v,„^ = ji^^M)dJ^ + ji^K.M) ^(^(-1)0^.), (3.21a) 



rP + v^Pp'^ = \^, + ,^ri _(,) , (3.21b) 



(1) , „,.7;(-i) _ \ „7; , ^ 

^n.«0f) = A,0, + --^^. (3.21c) 



When m = 1, the GKD of Kac-van Moerbeke equation is 



N 



Vt = t;(t;(-i) - t;(i)) + J^i^^'^'^^^j"^^^^'^ - vv^^^,4>f), (3.22a) 



i=l 



which has the Lax pair 

V'(^) + vtl)^-^^ = Xtfj, (3.23a) 

N ^2 _ _ AA - - 

+ E J^2^'^\''^M - ^^2^'^3h^^-'^- (3-23b) 

i=i -^ -^ 

The compatibility of the system (I3.23P under the condition (13.22b ) and (I3.22t i) 
will give (13.22b ) and (13.22b ). It is also an interesting question to find a Lax pair 
with /ij inside to give the system (I3.22p . This will be studied in future. When 
we take Aj = 0, the GKD of Kac-van Moerbeke equation (I3.22p will reduce to the 
Kupershmidt deformation version. When taking /Zj = 0, (I3.22p will lead to the 
Kac-van Moerbeke equation with self-consistent sources. 
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4 The generalized Kupershmidt deformation of 
Ablowitz-Ladik hierarchy 

Consider the following Ablowitz-Ladik discrete isospectral problem 

Z Q \ I _ I I I \T 



Ei; = U^, U=l^ ^^], ^ = (^,,^2)', (4.1) 

where Q = Q{n,t), R = R{n,t), z is the spectral parameter. Its adjoint problem 
reads 

E(-i)0 = 0f/, 0=(0,,02), (4.2) 

To derive the Ablowitz-Ladik hierarchy, we need first solve the discrete zero- curvature 
equation 

{Er)u - f/r = (4.3) 



Let 



we find that 



A B 
C -A 



(4.4) 



A^^h + B^^^R-Az-CQ = ^, (4.5a) 

A(i)Q + 5wi-5^ + Ag = 0, (4.5b) 

z 

C^i)^ - A^i? - Ai? - C- = 0, (4.5c) 

z 

C^^^Q - A(i) --BR + A- = Q. (4.5d) 

z z 

We can expand the above relations in power series of z and -, respectively [23]. 



4.1 r expanded in power series of I/2; 
Assume 

^ = Z^ r r-2*-i A --2i ■ (4.6) 

Then, the recursion relation (14.51) leads to 

4')-Ao = i?i = g(4') + Ao) Cf) = i?(4') + Ao), (4.7a) 

A« - A2, = gc2._i - i?i?£i = gci^li - i?i?2^+i, (4.7b) 

52m=Q(4l^+^2.) + 5£i, (4.7c) 

"2i+l 



4;|i = i?(A« + A2,) + C2._i, ^ = 1, 2, • ■ • (4.7d) 



where the initial value can be chosen as Aq = ^ 



2' 
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The Ablowitz-Ladik hierarchy can be written in the Hamiltonian form [23] 

-n. /, OU OU 

' I'm 

where Hm = —^^^ and the Hamiltonian operators are 

><''■■' - ( t t ) • ("^) 

Kn = QD-'Q^'\l - RQ) + (1 - RQ)Q^'^D~'EQ, 
Ku = E{1 - RQ) - QD-^R^^^E^'^\l - RQ) - (1 - RQ)Q^^^ D'^ER, 
K21 = -E'^~^\1 - RQ) - RD^^^^QE^^^\l - RQ) + {RQ - 1)R'^~^^D~^Q, 
K22 = RD^^RE{1 - RQ) + (1 - RQ)R^^^^D^^R. 

It is not difficult to find that ^ = ^20i, || = ^^i02- The GKD of Ablowitz-Ladik 
hierarchy is 

Q\ ^/-)f',y/,^-prtf^^^ (4.10a) 



We assume that 



Eipij = Zjipij + Qip2,j + fi,j, (4.11a) 

E2P2,j = Ri'l,j + -^2 J + f2,j, (4.11b) 

e(-^)0ij = Zj4)ij + i?02,i + gi,j, (4.11c) 

E(-l)02,i = Q4>1,J + -02,i + ^2,i. (4. lid) 



^i 



Simplifying (I4.10bp we get 



QD-'Q^^\l - RQ)^P2jhj + (1 - RQ)Q^'^D~^EQij2,jk,j 
+ Eil - RQ)i,ii4)2,j - QD-^R^^'^E^''\l - RQ)^ij4>2,j 

- (1 - RQ)Q'^^^D~^ER^^^,^2,, - z]{l - RQ)i>ij4)2,j = 0, (4.12a) 
E^-'\l - RQ)^2,jh, + RD-^QE^-^\l - RQ)^2,j4>ij 

+ (1 - RQ)R'^~^^D-'Qij2,j4>i,j - RD"^RE{1 - i?Q)^i,,02j 

- (1 - i?g)i?(^^)D-^i?^ij02j - ^1(1 - RQ)Mi,j = 0- (4-12b) 
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Substituting (14. lip into the above equation we find that 

z^{l - RQ)i4>2jfij - ^g^g) + (1 - /2g)g(^)[^S']^S!] + D-'Ei^,,g,j - 0i,,AJ] 

+ QD~'[{1 - RQ)i^2AQ^'^ ai]] - ^9i]i) - (1 - i?g)(^)02(^/i,. - ^jf2,j)] = 0, 

(4.13a) 

z,{l - RQ){i'2,m, - 4>[ff2f + (1 - RQ)R^~'\4>[7^fi:j'^ + D~\i^,,g,, - 4>iji,j)] 

+ RD-'[{1 - RQ)i^2AQ'^'^9^i.j - S^S) - (1 - RQY'^4>^2^{Rfi,, - z,h,)] = 0. 

(4.13b) 



First we assume that 



/2,i = -/ij, (4.14a) 



Zj 



92,j = -gi,j, (4.14b) 



then the above relations are simphfied and we can find that 

/i. = . ^\.. , (4.15a) 

91, = , ^' J , (4.15b) 

where /ij are some arbitrary constants. So the GKD of the Ablowitz-Ladik hierarchy 
reads 

EV^i,,- = z.V^ij + gV^2,i + T '^^^-. (4.16b) 

Eij2,j = R'4'i,j + -^2 J + . ^'^ . , (4.16c) 

Zj ^ (A ^j<t>2,j \ 



E(-^)0i, = z,4>i, + Rkj + . ,r , (4.16d) 

E^-^^2,j = Q4>i,j + -kj + ~ ^'^ 7 • (4.16e) 

W- 



Zj ' y (J, _ fi^^ 



Its Lax representation is 

E^jj = Uip, (4.17a) 

N 
^t = (Kn + J]X,-)^, (4.17b) 



J=l 



12 



.m, , / 



where Vm = {Tz'"')+ + , . 4 



X ' 



1 y2 _ y2 

Z Zj 









4 







2 Q 2 





J 



When we take fij = 0, the generahzed system (14.161) will reduce to the Ablowitz- 
Ladik hierarchy with self-consistent sources (corresponding to the generating matrix 
r in the power series of l/z). 

4.2 r expanded in power series of z 

Similarly, Assume F in ( 14.4p as 

^-[C -A -1^[C2.^,Z^^^' -A2.z^^ J' ^^-^^^ 

then we have the following recursion relations: 

4')-Ao = 0, fif) = -g(4')+Ao), C'i = -i?(4')+Ao), (4.19a) 

A^P - ^24 = gC'gli - RB2.-1 = QC2^+l - RB^]li, (4.19b) 

E« 1 = -g(4l) + A,.^ + 52,.i, (4.19c) 

C2i+, = -R{A'iJ + A2, + C^,, 2 = 1, 2, • ■ • . (4.19d) 

where the initial value is chosen as Aq = |- 

The Ablowitz-Ladik hierarchy can be written in the Hamiltonian form [23] 

-fL /, OM OU 

where ff^ = -%^5 -^^^^^ is given by (H^ . 

'^"^' - ( t; I:: ) ■ («i»' 

Kn = -QD-^QE{1 - RQ) - (1 - RQ)Q^^^D-^Q, 
K12 = E^^^\l - RQ) + QD'^RE^~^\l - RQ) + (1 - RQ)Q^-^^ D'^ R, 
K21 = -E{1 - RQ) + RD-^Q^^^E^^\l - RQ) + (1 - RQ)R^^^D'^EQ, 
K22 = -RD-^R^^\l - RQ) - (1 - RQ)R^^^D-^ER. 
So we have another type of the GKD of Ablowitz-Ladik hierarchy 

i^K^AL) _ 1 ^(AL)) f ^2,,01,, A Q_ ^4 22b) 
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k 


j ~ 


zQff, 




9i, 


j ~ 


zR-g], 


k. 


i ~ 




fij 


h. 


1 'Pi.j 

•1 "^ ZjRW 




i ~ 




/ij 


92,, 


V's 


1 ^iJ 



We assume that 

Ei^ij = Zji^ij + Qi}2j + fij, (4.23a) 

E^2,, = R^i,j + -^2 J + kj, (4-23b) 

E^-^^ij = zj^ij + R4>2,j + 9i,j, (4.23c) 

E(-1)02J = Q01,, + -02,i + ^2,i. (4.23d) 

Similarly with the case of F expanded in power series of 1/z, with some simplification 
we assume that 

(4.24a) 
(4.24b) 

then we can find that 

(4.25a) 

(4.25b) 

(4.25c) 

where Hj are some arbitrary constants. So another type of the GKD of the Ablowitz- 
Ladik hierarchy is 

EtPij = z,ij,j + Qi/J2j + _ ^'^'^^^ . , (4.26b) 

E^2,j = R^i,j + -^2,j + ^^^^, (4.26c) 

^ ^^'3 ^ z,RW 

^("1)01, = z,0i,, + R^2, + _ ^'''^. . , (4.26d) 

E(~'^2,j = g0i, + ^02, + _ ^%^^ • (4.26e) 

Its Lax representation is 

E^jj = Uip, (4.27a) 

N 

V^t = (K. + J]^,)^, (4.27b) 

i=i 
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2m\ I / —■^2m. 



where Ki = (r^; /- . , g q 



X 



' z^ - z] 



ZZj1p2,j(Pij Z 1P2,J<P2,J 



4 







2 o 2 
2: — o^ 








When we choose ^j = 0, the generahzed system ( J4.26P will reduce to the 
Ablowitz-Ladik hierarchy with self-consistent sources (corresponding to the gen- 
erating matrix F in the power series of z). So we get two types of generalized Ku- 
pershmidt deformation of the Ablowitz-Ladik hierarchy according to the two types 
of Hamiltonian operators (14.91) and (I4.2ip . 



5 Conclusion 

In this paper, we construct some new integrable discrete systems and their Lax 
representations by making use of the generalized Kupershmidt deformation of bi- 
Hamiltonian systems. The generalized Kupershmidt deformation for the Toda hier- 
archy, the Kac-van Moerbeke hierarchy and the Ablowitz-Ladik hierarchy are stud- 
ied. We find that these new discrete systems have some relations to the equation 
with self-consistent sources [H [3 [TJ [20] . It is shown that the method can be used to 
construct new discrete integrable systems in (1 + l)-dimensional case. The Kuper- 
shmidt deformation for higher dimensional case still needs to be studied. 
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